The structural phase transition and loss of magnetic moments in NpC^: 
ab initio approach to the crystal and mean field 
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We discuss the triple-g* x structures for the ordered cubic phase of Np02, which are Pn3m and Pa3. 
A special care should be taken to discriminate between these two cases. We analyze the relevant 
structure-factor amplitudes and the effect of domains on resonant X-ray scattering experiments. We 
formulate the technique of multipole expansion of the Coulomb interaction and use it to study the 
crystal electric field and mean field for a number of neptunium many electron configurations (5/ 3 , 
7s5/ 3 , 7p5/ 3 and 6d5/ 3 ) on ab initio level. We have found that the crystal field is rather small 
(50-150 K). The direct quadrupole-quadrupole interaction between neighboring neptunium sites is 
weak and can not drive the structural phase transition at 25 K. We have introduced an effective 
(enhanced) quadrupole interaction and considered the interplay between it and the crystal electric 
field. The influence of both interactions on the transition temperature has been investigated in 
detail for the 5/ 3 configuration. We discuss the importance of the intrasite multipole interaction 
between three localized electrons (5/ 3 ) and a valence electron on a neptunium site. We show that 
this interaction combined with the symmetry lowering at 25 K may be responsible for the loss of 
the magnetic moments in the ordered phase of Np02- 

PACS numbers: 64.70.Kb, 71.70.-d, 71.70.Ch, 75.10.Dg 



I. INTRODUCTION 

Since the discovery of the low temperature phase tran- 
sition in NpC>2 at T c — 25.5 K almost fifty years agoji 
the magnetic properties of this compound represent a 
challenge for theoretical interpretation!^ NpC>2 and its 
neighbor in the actinide series, UC^r^ are both crystal- 
lized in the cubic fluorite (CaF 2 ) structure in the high 
temperature phase. The temperature dependence of the 
magnetic susceptibilities of these compounds is sugges- 
tive for a Neel transition to an antiferromagnetic state4£ 
However, while it was found that UO2 orders antiferro- 
magnetically (Tn — 30.8 K) with 1.74/xb/U atom^ no 
magnetic ordering^ was confirmed for NpC>2 (see for a 
review Ref. @) • As follows from Mossbauer experiments^ 
the upper limit on the magnetic moments of Np in the 
ordered phase is only 0.01/zs/Np. This implies almost 
a complete quenching of the 3.0/is magnetic moment of 
Np, which it exhibits in the paramagnetic high temper- 
ature phase. The demagnetization of Np magnetic mo- 
ments seems to be incompatible with the widely accepted 
viewpoint that the Np cations are in the tetravalent va- 
lence state, Np 4+ , with three localized 5f electrons left 
at each neptunium site. With the 5/ 3 configuration, Np 
becomes a Kramers ion having a magnetic doublet (Tq, 
Tr) or quadruplet (Fg) ground state. Because of this, a 
divergent susceptibility at T — > will be present in any 
model of the phase transitionAS 

To describe the disappearance of magnetic moments of 
Np at T < T c Santini and Amoretti put forward an idea 
of a magnetic octupole order parameter,^ which on one 
hand, is not invariant under time reversal symmetry and, 
on the other hand, is different from the magnetic dipolar 
order parameter which brings about the ordinary mag- 
netic ordering. The idea of breaking down of Kramers' 



degeneracy by the magnetic octupole seems to be the 
only solution in the framework of the 5/ 3 modeli^ 

An alternative explanation was offered by Friedt et 
al&. They suggested that the magnetic anomaly in NpC>2 
could be explained by a splitting of the ground state 
T$ quadruplet of cubic symmetry into two doublets in 
the crystal field of lower symmetry of the ordered phase. 
The splitting was ascribed to an internal distortion of the 
oxygen sublattice at T c . This viewpoint was inspired by 
the experimental observation that the phase transition 
in UO2 is accompanied by an internal distortion of the 
oxygen cube that surrounds the U cation^ while the ex- 
ternal cubic structure of uranium dioxides survives the 
transition. However, the mechanism mainly weakens the 
magnetic response from the Np sites at 20 K, while the 
T — > divergence of magnetic susceptibility still persists 
in the model in contradistinction with experiment On 
the other hand, no evidence for an internal or external 
crystallographic distortion in NpC>2 has been found by 
synchrotron experiments, 11 Thus, the phase transition 
in NpC<2 appears to be isostructural like the 7 — a phase 
transition in elemental cerium 1 ^ 3 . (T c ~ 100 K) or the 
isostructural expansion in YbInCu4 at T = 42 K. 14 It is 
worth mentioning that in all these compounds the phase 
transition is accompanied by a loss of magnetic moments 
in the ordered phase. 

However, very recent resonant X-ray scattering (RXS) 
experiments at the Np Miy and My edges in NpC-2 
indicated an unexpected result: the phase transition is 
not isostructural^ In the low temperature phase a long 
range order of Np electric quadrupoles was revealed by 
the growth of superlattice Bragg peaksii 1 * 1 ^ The space 
symmetry of the ordered phase was identified as Pn3mm^ 
The symmetry lowering is a special one. The centers of 
mass positions of neptunium and oxygen remain in the 
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cubic CaF2 structure as in the high temperature phase. 
However, electronic quadrupoles of the Np sites have four 
different orientations, which allows us to distinguish four 
different sublattices of Np cations^ 



This experimental finding gives rise to a question of 
correlations between structural and magnetic properties 
in solids, which we have investigated theoretically in our 
model of the 7 — a phase transition in Cem^H^ Unlike 
Np02, pristine cerium is a metal and the phase transition 
there at normal pressure is accompanied by a huge vol- 
ume change^2*i£ The volume anomaly in NpC>2 (contrac- 
tion) is only 0. 01 8% However, the disappearance of the 
magnetic moments and the "isostructural" character of 
the phase transition makes them similar. In this respect 
it is interesting to notice that in our model for Ce we 
have predicted the Pa3 space symmetry for the ordered 
a phase 16 which is very close to the Pn3m structure re- 
ported for NpC>2 in Ref. The active irreducible rep- 
resentation also belongs to the X point of the Brillouin 
zone (BZ). While the existence of superstructure reflec- 
tions for a-Ce remains an open question which has to be 
investigated experimentally^iLi^ here we want to ap- 
ply our theoretical concepts for the study of the crystal- 
and mean-field in neptunium dioxide. 



We will use the technique of multipole expansion of 
the Coulomb interactioniiSiiLiSiiii The multipole expan- 
sion represents a unified description of the crystal field 
effects and atomic term splitting. The intrasite Coulomb 
repulsion which is responsible for Hund's rules, the spin- 
orbit coupling and the crystal field effects are included on 
equal footing. (For a single site the technique is equiv- 
alent to the classical description of atomic termsi^) In 
order to calculate the corresponding electronic spectra 
we use many Slater determinants, which indicates that 
our scheme is a genuine many electron approachi&iii cor- 
responding to the configuration interaction (CI). (Notice 
that ordinarily used single-determinant Hartree-Fock ap- 
proach is not sufficient for these purposes). 



The paper is organized as follows. In Sec. II we discuss 
the triple-g* structures (PnSm or Pa3) of the ordered 
phases and derive the selection rules for resonant X-ray 
scattering (RXS) experiments. In Sec. Ill we describe our 
method for treating local electron configurations at nep- 
tunium sites. Next (Sec. IV) we examine thoroughly the 
case of three localized 5/ electrons. In Sec. V we general- 
ize the 5/ 3 model by including the multipole interactions 
with a valence electron which is present instantaneously 
on a Np site. We show that in the framework of this 
extended model the disappearance of the magnetic mo- 
ments can be understood and ascribed to the trigonal 
symmetry lowering at 25 K. Our conclusions are summa- 
rized in Sec. VI. 



II. STRUCTURAL PHASE TRANSITIONS AND 
RXS EXPERIMENTS 

In this section we discuss the allowed space groups of 
the low temperature phase of Np02- Our considerations 
here are based on a general group theoretical approach 
which is independent of the model assumptions. 

In the disordered phase (T > 25 K) the crystal struc- 
ture is cubic {FmZm) and the site group is Oh- The 
electron density of neptunium is given by the spheri- 
cally symmetric component (V °) and the cubic harmon- 
ics Ki{Q) with I — 4,6. There is no contribution to 
the density from the quadrupole functions Y^_ 2 i where 
t = 0, (m, c) or (m, s), m = 1,2. (Here and below we 
work with the real spherical harmonics with the phase 
definition of Ref. iJ) At T c = 25 K the transition to 
a new phase sets in. The new phase was characterized 
as a triple-<f antiferro- quadrupolar ordering of T 2g (r5 s ) 
electric quadrupoles at the Np sites^ 

In real space the ordering is characterized by four 
different sublattices of the simple cubic structure. We 
label these sublattices which contain the sites (0,0,0) 
(o/2)(0,l,l). (o/2)(l,0,l) and (o/2)(l,l,0) by {n p }, p = 
1 — 4, respectively. The most significant feature of the 
ordered phase is the existence of only one three-fold axis 
of symmetry C3 at each Np site which is also a cube di- 
agonal. The only quadrupole function compatible with 
the symmetry lowering is Y 2 (f2') in the coordinate sys- 
tem where the z'-axis coincides with one of the three-fold 
axes (cube diagonals): [111], [—1,-1,1], [1,-1,-1], and 
[—1,1,-1]. Consequently, there are four such functions 
which are given by 



5„(n) = 








(2.1a) 


S b (Cl) = 




- y^m 




(2.1b) 


S C (Q) = 






Y 2 2s m, 


(2.1c) 


S d (Sl) = 




+ Y 2 lc (Q) 


-Y 2 2s m. 


(2.1d) 



The spherical harmonics Y 2 ls , Y 2 lc , Y 2 2s belong to a three- 
dimensional irreducible representation T 2g of Oh- They 
are proportional to the Cartesian components yz, zx, and 
xy. 

In Ref. [jja the space group of the ordered phase was 
identified as Pn3m (No. 224, Ref. l22j) . The correspond- 
ing ordering of S a — Sd functions is shown in Fig. ^ left 
panel. However, a close examination reveals that there 
is another possibility which was not considered by the 
authors of Ref. This is the Pal structure (No. 205 of 
Ref. 113) depicted in the right panel of Fig. [TJ The Pa3 
structure also realizes a triple-cf quadrupolar order, the 
site symmetry of the Np sites is = C3 X i. 

The problem of finding the appropriate space group 
can be simplified to the task of distributing four orien- 
tational functions S a — Sd among four distinct sublat- 
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Pn3m Pa3 

FIG. 1: Pn3m and Pa3 structures for quadrupoles of Np02. 
Both structures have the common Se (or 3) local site symme- 
try. They differ in the way the threefold axes and quadrupoles 
are distributed over 4 sublattices {n p }, p = 1 — 4, (see text for 
details). The Pn3m structure possesses three mirror planes 
one of which (m) is shown explicitly. 

tices {Hp}, p = 1—4. Keeping the three-fold rotation 
symmetry about the [111] axis and S a for the first sub- 
lattice, we distinguish only three different choices here. 
The first is when (i) Sb corresponds to {714}, S c — {n 2 }, 
and Sd~- {n 3 }. This gives the Pn3m structure, Fig.Q] left 
panel. Another assignment is (ii) Sb — {^2}, S c — {^3}, 
and Sd~- {^4}. This scheme corresponds to the Pa3 struc- 
ture as shown in Fig. ^ right panel. The third choice is 
(Hi) Sb — {n 3 }, S c — {^4}, and Sd — {n 2 }. This is simply 
another variant (or domain) of the PaZ structure. 

From the mathematical point of view we are analyz- 
ing the symmetry lowering and the condensation schemes 
at the X point of the BZ of the fee lattice, which in- 
volves the density components of the T 2g symmetry. In 
both cases (Pn3m and Pa3) the triple-<f mode which 
drives the structural phase transition belongs to the X 
point of BZ *q x , and involves the three arms of the 
latter, qf = (27r/a)(l, 0, 0), q x = (27r/a)(0, 1,0), and 
q x = (27r/a)(0,0, 1). To describe the three quadrupolar 
components S\ = Y 2 ls , S 2 = Y 2 lc , S 3 — Y 2 2s at site n we 
introduce the functions Si(n) (i = 1 — 3) and consider 
their Fourier transforms, 

n 

where the position vector X runs over N sites of the face 
centered cubic (fee) neptunium lattice. The little group 
of *q x is Dm (A/mmm). At q x the functions S 2 and S3 
belong to the E g small representation, while the function 
5*1 belong to the B 2g small representation. Consequently, 
we distinguish two different irreducible representations of 
the space group FmSm. The first one, X~£ in notations 
of Stokes and HatchjS 3 . comprises six orientational func- 
tions, with two functions from every arm of *q x . The 
functions are S 2 (q x ), S 3 (q x ); S^qf), S 3 (q x ); S^qf), 
and S 2 (q x ). The second irreducible representation, X£, 



has only three components, with one function from each 
arm of *q x , namely Si(q x ), S 2 (q x ) and S 3 (q x ). For the 
symmetry lowering to the Pn3m structure we obtain the 
following condensation scheme: 

Fm3m : X^ 

((Siffi)) = (S 2 (q x )) = (S 3 (q x )) = Pl VN £ 0) 
-> Pn3m (Z = 4). (2.3) 

Here (...) stand for a quantum and thermal average and 
pi is the order parameter amplitude. The X± active 
representation corresponds to the irreducible star {fcio} 
and the loaded one-dimensional representation f 7 in Ko- 
valev's notation^ There are four domains for the Pn3m 
structure. 

The condensation scheme for the Pa3 structure is given 

by 

Fm3m : X^ 

i(Ss(q x )) = - (5 2 (gf )) = P 2^N ? 0, 

(S 2 (q' X )) = (S 3 (q X )) = (S 1 (q X ))=Q) 

-> Pai (Z = 4) domain I. (2.4) 

Here again p 2 is the order parameter amplitude. This 
condensation scheme corresponds to a domain shown in 
Fig. 2] and there are eight possible domains of Pa3. In 
particular, there is a domain (iii) which we have already 
considered before (<S Q — {ni}, Sb — {n 3 }, S c — {Hi}, and 
Sd — {n 2 }) when we analyzed the variants of distribution 
of S a — Sd over four sublattices. This domain is obtained 
as a result of the following symmetry breaking 

Fm3m : Xg 

((S 3 (qX)) = (Si(q X )) = (S 2 (q x ))=0, 

= (Ss(tf)) = (Si(q x )) = P2 ^N ± 0) 
-> Pal [Z = 4) domain II. (2.5) 

Previously the symmetry change Fm3m — ► PaS has been 
considered in a number of publications, Refs. H^IT^sl 
[2(J The X§ active representation is defined by the star 
{kio} and the loaded two-dimensional representation f 9 
in Kovalev's notation^ 

Although the two structures look similar there is 
one very important difference between them. The Pa2> 
ordering corresponds to an effective attraction of the 
quadrupoles, while the Pn3m ordering results in a re- 
pulsive interaction between them, see Appendix A for de- 
tails. This is the main reason why the PaZ structure (but 
not Pn3m) is found in many molecular solids such as H2, 
N 2 £L NaC-2^ and C 60 . 26 i 28 i 29 i 30 i 31 Here we speak about 
the direct bilinear electronic quadrupole-quadrupole in- 
teraction. However, in case of NpG<2 there are indirect 
(superexchangei 2 ") interactions via oxygen which may 
lead to an effective attraction and therefore, the Pn3m 
structure can not be ruled out. Nevertheless, the Pa3 
space group is a good candidate for the symmetry of the 
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ordered phase of NpC>2- Provided that the superexchange 
interaction remains the same for two structures, the Pa3 
symmetry becomes preferable. It is interesting to notice 
that at first the symmetry of the ordered phase of the 
pristine Cgo molecular crystal was identified as Pn2>?- 
Only later it was corrected and proved to be Pa3^^^ 
Now we will study how the Pn3m and Pa3 struc- 
tures manifest themselves in resonant X-ray scatter- 



ing experiments 



15,33 



The tensor of scattering on the 



quadrupolar density of a neptunium atom in a cubic lat- 
tice is given by 34 i 35 i 36 



pi p2' 

fn = f | Pi p 3 
K P2 P3 



(2.6) 



where pi = ±1. We have the following 
correspondenc e - 35 ' 36 between the functions S n , Eq. I|2.1at 
d), and the scattering tensors: the function S a 
corresponds to f a , where pi = p2 = pa = 1, the function 
St to fi, with pi = 1, p2 = p3 = — 1; the function S c 
to f c with pi = p 2 = -1, P3 = 1; and S d to f d with 
p± = P3 = — 1, P2 = 1- The scattering tensor structure 
amplitude is found as 



N 



F(h,k,l) = J2fn 



(2.7) 



where K = (2Tr/a)(h,k,l) stands for the vectors of the 
reciprocal lattice and a is the cubic lattice constant. In 
performing the summation in 1|2.7|) we distinguish four 
contributions from four sublattices {n p }, p = 1 — 4. We 
are mainly interested in superstructure Bragg reflections. 
We obtain then that in general 

T(h,k,l)=FM, (2.8) 

where F = fN and the matrix M. is either A, B, or C: 



A = 




B = 




C = 




(2.9) 

In Table [I] we quote which of the matrices (A, B,C) oc- 
curs for every particular case of (h, fc, I) for the Pnim 
and Pa3 structures. We recall that the conditions^ for 
the isotropic scattering on the fee lattice of Np (from the 
spherically symmetric densities) are h + fc, fc + I = 2n, 
that corresponds to two first lines of Table Q] 

Finally, we consider the polarization dependencies for 
resonant X-ray scattering experiments. We assume that 
a crystal has a flat (001) surface and the azimuthal angle 
-0 is counted from the a:— axis defined by e x — (1,0,0)^ 
The incident beam direction is given by (—cos 9 cos?/', 
— cos 6 sintp, — sin 6), the scattered beam direction by 
(— cosG cosV>, —cos 8 smip, sin0). We introduce the 



TABLE I: Tensor T(h, k, I) = FM for the superstructure 
Bragg reflections on quadrupolar densities of Np. A4 = A, B 
or C, Eq. 112. 91 . I (ii) and II (Hi) refer to two domains of the 
Pa3 structure described in the text, is the zero matrix. 



Bragg reflection 
h,k,l : 




Pn3m 

(0 


I(ii) 


roo 

II (Hi) 


h,k,l = 2n 













h, k,l = 2n+ 1 













h = 2n + l, k,l = 


2n 


A 


c 


B 


k = 2n + l, h, I = 


2n 


B 


A 


C 


l = 2n + l, h,k = 


2n 


C 


B 


A 


h,k = 2n + l, 1 = 


2n 


C 


B 


A 


h, I = 2n + 1, k = 


2n 


B 


A 


C 


k,l = 2n + l, h = 


2 n 


A 


C 


B 



TABLE II: 


The structure factor amplitudes F p ^ p /(h, k,l)/F 


for four polarization channels of RXS and the matrices A, B 


and C. 










A 


B 


C 


a — > a 








- sin(2V0 


a^iYf 


cos O COS ij) 


— cos sin ip 


cos(2i/>) sin 


iTi — > a 


cos Q cos ij) 


— cos Q sin ip 


— cos(2V0 sin 


TVi -> IV f 








sin(2V>) sin 2 



standard polarization vectors e*i (i = 1 — 3) parallel (tTj, 
TTf) or perpendicular (a) to the scattered plane, 

a = (— sin?/', cosV', 0), (2.10a) 
7?,; = (— sin6 cos?/;, — sin sin?/', cos0), (2.10b) 
iff = (sin0cos^, sin0sin?/', cosO). (2.10c) 

A structure factor amplitude F(h, fc, I) is found a^i^iS 6 . 



F p ^ p , (h, k,l) = F e p , T M e p 



(2.11) 

where T stands for transpose and p — > p' denotes one of 
the four polarization channels, a — > a, a — > ttj, ii{ — > a, 
and TTi — > nf. The structure factor amplitudes for the 
three cases (A4 = A,B,C) are quoted in Table [H] The 
corresponding intensities of the superstructure Bragg re- 
flections are found as I p ^, p i(h,k,l) = \F p ^ p i (h, fc, l)\ 2 . 

From TablesmandlTTlone can easily obtain all necessary 
dependencies for intensities of different polarizations. For 
example, the intensity of the (003) reflection of Pn3m is 
exactly the same as the intensity of the same polarization 
of the (300) reflection of domain I (ii) of Pa3, and the 
intensity of the (030) reflection of domain II (Hi) of Pa3. 
Furthermore, the (013) reflection of domain I (ii), and 
the (103) reflection of domain II (Hi) of Pa3 also have 
the same intensity, and so on. Since both symmetries are 
very similar, we believe that special care should be taken 
to discriminate between the Pn3m and Pa3 structures. 

Tables HI and ITU are also useful in considering the con- 
tributions from the domains of the same group. We start 
with the Pn3m structure. The four domains differ by 
the orientation of the y 2 °(f2') quadrupolar function at 
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n = 0. This function may be chosen to align along 
four main cube diagonals, which are connected with each 
other through rotations by ir/2 about the z— axis. Hence, 
all four domains are obtained from the first, (i), by apply- 
ing three consecutive rotations by 7r/2 about the z— axis. 
The effect of the domains in RXS experiments is equiva- 
lent to a superposition of four structure amplitudes, Ta- 
ble^ with angles -0, + 7r/2, ip + n, and + 371-/2. (Here 
we suggest that the domains produce the coherent scat- 
tering.) We label these domains by indices dl — dA, and 
introduce their populations, Pdi — Pdi- The population 
of a domain d is defined as Pd = Nd/N, where Nd is the 
total number of neptunium atoms in the domain. Then 
we obtain from Eq. (|2.11|l that the structure amplitude 
for the a — ► a and a — > 7T/ channels of the (003) reflection 
is given by 

^-^(003) = -PP sin(20), (2.12a) 
F ff ^ 7r/ (003) = PP cos(20) sin 6, (2.12b) 

where 

P = {Pdi ~ Pd2 + Pdz ~ Pdi) ■ (2.12c) 

If Pdi + P d3 = Pd2 + Pdi, P = and the (003) 
Bragg reflection is suppressed. Otherwise, the ip— and 
9— dependencies in (|2.12al b~) are exactly the same as for 
a sing le domain, Tables EHO For the (300) reflection of 
Pn3m one obtains that 

F^„ f (300) = PP' cos9cos(^ + o ), (2.13a) 

where 

= arccos([P d i - P d3 ]/P'), (2.13b) 

P' = y/[P d i - Pds] 2 + [Pd2 ~ Pdi] 2 - (2.13c) 

Comparing this result with the expression for a single do- 
main, we observe that the main effect is the phase shift 
ipo given by Eq. I|2.13b(l . The condition for suppression 
of the reflection is Pdi = Pd3 and Pdi — Pdi- The po- 
larization dependencies of the other reflections of PnZm 
can be figured out analogously. 

There are eight domains of Pa3 structure. Earlier we 
have already considered the two basic variants of Pa3: 
I (ii), and II (Hi). The others are obtained by rotating 
these two variants by the angles 7r/2, 7t, and 37r/2 about 
the z— axis. Applying the rotations to I and II, we arrive 
at two series. We label the population of the series of 
domains by the indices P\ x — Pj 4 , and Pjl — PjJ , respec- 
tively. The two series result in two distinct scattering 
matrices, for example, for (003) they are B and A, for 
(300) C and B, etc, see Table [I] The domain pattern 
produces a superposition of 8 terms (amplitudes). Each 
term corresponds to one of the two matrices and to one 
of the four azimuthal angles: ip, ip + tt/2, ip + n, and 
tp + 37r/2. The relevant expressions are obtained in the 
same way which we have used to derive Eqs. (|2.12al -c') 



and l|2.13al -c'). For example, 

^^.(003) = 0, (2.14a) 

^^(003) = F cos 9(-P' / sin(-0 + i/jq) 

+P' 11 cosOj + rt/)), (2.14b) 

The latter expression can be transformed to a shifted 
sin— or cos— like functions, i.e. P cr ^ 7r/ (003) ~ 
cos 9 sin(^ +ip'o), where ifj' is a phase shift depending on 
the domain pattern of Pa3. P' 1 , P' 11 , and ipo, ^o 1 are 
given by Eq. (|2.13c|l and (|2.13bf) . where P d i (i = 1 — 4) 
are replaced by PL or PU . For the (300) reflection of 
Pa3 we find 

F(j — (300) F P^ sin (2^) , (2.15a) 
F a ^ nf (300) = FP 1 008(2?/-) sin 9 

-FP' H sin(ij + V'o / )cos9. (2.15b) 

Here again P 1 is given by 12.12c|) for PL, i = 1 — 4. 
As before, the superstructure Bragg reflections are sup- 
pressed in the case of equal population of 8 domains in 
the crystal. 

III. CONFIGURATION INTERACTION 
CALCULATION 

In this section we describe in detail how we construct 
the basis consisting of many determinantal wave func- 
tions for a many electron system and how we calculate 
the relevant matrix elements. First, we perform a band 
structure calculation in order to determine the conduc- 
tion electron charge distribution in the muffin-tin (MT) 
sphere around Np. Subsequently, we apply the config- 
uration interaction method to treat the many electron 
system consisting of the localized 5/ 3 electrons and the 
conduction electron in the MT sphere. 

A. Electron band structure calculation 

We have started by performing an electron band struc- 
ture calculation of Np02 using our linear augmented 
plane wave (LAPW) code^i The calculation has been 
done assuming the muffin-tin (MT) shape of the one- 
electron potential and the Barth-Hedin expression for 
exchange^ which is a variant of the local density ap- 
proximation (LDA). The equal MT radii R^ P T =2.2206 
and R^jrp =2.2206 in atomic units (a.u.) were chosen for 
Np and O, with the cubic lattice constant a = 10.2567 
a.u. (or 5.4276 A)pii The MT potential and density of 
Np and O have been obtained self-consistently using a 
LAPW basis of ~300 functions on a 20-point mesh of 
the irreducible part of the Brillouin zone. The three 5/ 
electrons of Np were treated as core states, which adjust 
self-consistently to the conduction electron density. 

As a result of the calculation we obtained that Np02 
is an insulator, with the energy gap AE = 0.789 eV. 
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TABLE III: Angular-momentum-decomposed (partial) elec- 
tronic charges Qf and total charges Q A inside neptunium 
and oxygen MT spheres and in the interstitial region (LAPW 
calculations, see Ref. l39l for details and definitions); Q Np = 
-(2Q°+Q l ). 



A 


Np 





interstices/per unit 


Qi 


0.041e 


0.027e 






0.401e 


4.313e 




Qd 


0.261e 


0.046e 




Qf 


2.185e 


O.Olle 




# 


+4.108|e| 


-0.284je| 


-3.540|e| 



The width of the occupied electron band below the Fermi 
level is Et — Eb = 5.953 eV. The spin-orbit splitting be- 
tween 5/7/ 2 and 5/5/2 one-electron states is A/ =0.983 
eV. The main goal of the calculation however was the 
electron charge density distribution inside the neptunium 
MT-sphere. The calculated partial charges of different 
angular symmetry (Z = — 3) are quoted in Table II I II 
An important result is that on average there is approx- 
imately one conduction electron present inside the Np 
MT-sphere. Therefore, the localized 5/ 3 configuration of 
Np can not be considered separately from this valence 
electron, which can be of 7s, 7p or Qd type. The in- 
stantaneous configuration at the neptunium site becomes 
7s5/ 3 , 7p5/ 3 or 6d5/ 3 . 

In all cases this additional electron experiences a strong 
Coulomb repulsion with the three localized partners. 
This interaction is not fully accounted for by the band 
structure calculations because it requires a multideter- 
minant treatment or configuration interaction (CI). 18 
Therefore, we have to follow a different route and be- 
low we study the electron spectrum using the multipole 
expansion of electronic densities. 



B. Many electron basis states 

Our method of multipole expansion of the Coulomb 
interaction has been used beforeiiSiiLiLiS Here we for- 
mula te it i n detail for the s/ 3 configuration following 
Refs. Il8ll9i In Sec. IV we deal with the 5/ 3 configura- 
tion which is easily obtained from s/ 3 by omitting one s 
electron. In Sec. V we consider 7s5/ 3 , 7p5/ 3 and 6d5f 3 
configurations. For the latter two cases we will describe 
the important differences with s/ 3 . 

We start by considering a face centered cubic (fee) crys- 
tal of N Np atoms. Each atomic site possesses one 7s and 
three 5/ electrons. The position vector of an electron 
near a crystal lattice site n is given by 

R(n) = X(n) + f(n). (3.1) 

Here X(n) is the lattice vector which specifies the centers 
of the atoms (Np-nuclei) on a rigid fee lattice. The radius 
vector f(n) is given in polar coordinates by (r(n), f2(n)), 



where r is the length and f2 = (0, <fi) stands for the polar 
angles. We label the basis ket-vectors at the lattice site n 
by a single index I or, alternatively, by four one electron 
indices i£, 1$; i s )'- 

= (3.2) 

The index = (m*, sj) stands for the orbital (m? = 
1 — 7) and spin projection (s z = ±1/2) quantum numbers 
of one / electron. Therefore, there are 14 states which 
we label by i* = 1 — 14. Two states of the s electron 
are labeled by i s — 1,2. The many electron basis wave 
functions are 

3 

xl[{?t\i f t)ri-m s )ii, (3.3) 
t=l 

where a stands for a permutation of four electrons, the 
factor (— l) p takes into account the parity of the permu- 
tation, N a is the number of the permutations, and 

(r\i f )n=n f (r(n))(h\if), (3.4a) 
(f'\i s )n = Tl s (r'(n)){n , \i s ). (3.4b) 

Here IZf and 1Z 8 are the radial components of the 5/ 
and 7s electrons, respectively; n stands for fl(n). The 5/ 
spin-orbitals can be written as 

(h\if) = {h\m f )u s (s z (f)), (3.5) 

Here u s is the spin function (s = ±). The /—orbital 
parts, (n|m/), are expressed in terms of spherical har- 
monics Y{ a (yi) = (n\l, m). We find it convenient to work 
with real spherical harmonics 21 Y ; T , where r = 0, (m,c) 
or (m, s). 

The order of indices in (|3.2II is important. For example, 
as follows from the dynamical equivalence of the electrons 
the state i{ , i s , i^) can be reduced to i s ) by 

permuting the third and the fourth electrons, i.e. 

|i{,4,i s ,z^) = -|z{,4,z^ ;l s ), (3.6) 

and so on. To describe the same quantum state we will 
use the basis vectors (|3.2|) and apply the permutation 
law (|3.6|l when needed. Alternatively, one can use the 
corresponding Slater determinants for the four electron 
wave functions, Eq. (|3.3|) . However, the permutation 
relations of the type of Eq. (|3.6I) are more efficient for 
our purposes. Excluding equivalent states, we find only 
(14 ■ 13 • 12/3!) x 2 = 728 independent functions, or de- 
terminants, for 7s5/ 3 . (These are 364, 2184 and 3640 for 
5/ 3 , 7p5/ 3 , 6d5/ 3 , respectively.) Notice, that every basis 
wave function is in fact a Slater determinant, Eq. (|3.3() . 

C. Multipole repulsion between electrons 

Now we take into account the Coulomb intrasite and 
intersite repulsion by expanding the interaction in multi- 
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pole series. As was discussed in Ref.0, these interactions 
are treated exactly in the chosen quantum space (7s5/ 3 ). 

The Coulomb interaction between two electrons at 
sites n and ft' is given by 



V(R(n),R'(n')) = 



1 



\R(n) - R'{n')\ 



(3.7) 



The multipole expansion in terms of site symmetry 
adapted functions (SAF's)2i is 



V(R(n),R'(H')) =J2 v AA>{fi,fi'; r,r') S A (h) S A ,{ti), 

(3.8) 



AA' 



where 

vaa' (n,n r ; r,r') 



dn(n) / dn'(n) 



Shin) S A >(n') 



\R(n) 



R'(ft')\ 
(3.9) 



The SAF's are linear combinations of spherical harmon- 
ics and transform as irreducible representations of a site 
point group, Ref. |2l[ The index A stands for (1,t), with 
t = (r, /i, k). Here I accounts for the angular dependence 
of the multipolar expansion, T denotes an irreducible rep- 
resentation (in the present case the group is Oh), p labels 
the representations that occur more than once and k de- 
notes the rows of a given representation. 

The intrasite case corresponds to n = nf. The in- 
teraction function VAA'(r, r') = VAA'{n = n' \ r,r') then 
becomes particular simple, 



v\A'(r,r') = 



-In 



21 + 1 



'AA' 



(3.10) 



where r> = max(r,r'), r < = min(r,r') and <5aA' = 
S TT 'Su'. The last expression is also site independent. 

There is no simple analytical expression for the inter- 
site case, n ^ n'4£ The intersite multipole interactions 
are anisotropic and for practical purposes it is important 
to use the following dependence^ 



v A a> (n,n'; r,r') 



{r) l {r') V 



\X(n) - X(n')\ l + l '+ 1 



(3.11) 



D. Intrasite matrix elements 

For the Coulomb interaction between four electrons on 
a same site n we have a sum of six two-body terms, 
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^ (4) = ^E 2J r(r " r «>- 



(3.12) 



where each term is given by the multipole expansion 
(|3.8|) . In order to calculate the matrix elements of 
(i(,i2,il; i s \V^\j{ij{,j{]j s ), one has to classify the 



electronic transitions. Following Ref. [l9| where the en- 
ergy terms of molecular ions C^, to = 2 — 5, were 
calculated, we consider four possibilities for the fourth 
s— electron: (1) i s — > j s , (2) i s — > jg, (3) i s — > j[ , and 
(4) i s — > j[, which we label by the index a± = 1 — 4. 
The et4 = 2 and 04 = 4 transitions involve odd number 
of transpositions among j( , j| , j(; j s , and the parity is 
P(a4 — 2) = P(d4 = 4) = —1. For two other transitions 
the number of transpositions is even and P(a4 = 1) = 
P(a4 = 3) = 1. After this we are left with only three 
j-states, which we label as j[, j' 2: and j 3 . For the next 
electron, ig, we can consider three possibilities (ig — > 

*3 ~~ * J2' *3 ~~ * Ji) which we label by the index 03 = 1 — 3. 
In this way we continue until we exhaust all four elec- 
trons. As a result, each subcase (or electron transition) 
is classified by the three index label a = (04, 03, (12), and 
its parity is P(a) — P(a4)P{a^)P{a2)- Mathematically, 
we reduce a permutation of four electrons to a product of 
transpositions. The matrix element (I\V^\J) is found 



(I\V {4) \J) =J2 p ( a )( I \V {4) \J) {a) , (3.13) 



where £ Q = Et t =i E! 3 =i £a 3 =i> and 

l,T 



><S(4,ja2)S(i{,jal) + P-i 



f 



(3.14) 



Here p.i. stands for the other pair Coulomb interactions, 
Eq. (|3.12|l . (The explicitly written term in Eq. (|3.14|) 
corresponds to the interaction between fourth and third 
electron.) The elements ca(?j) = Q, T (y) are defined by 



c A {ij) = dQ.(i\h)S A (n)(n\j). 



(3.15) 



For the 7s5/ 3 configuration there are three types of these 
coefficients. For the s — s transition it is only one in- 
tegral (s|Yo |s) = I/V^tt, which is not zero. For the 
/ — / transitions and real spherical harmonics^ the co- 
efficients c^ r (if j?) were tabulated in Ref. [l|| Finally, 
there are f — s and s — f transitions which require the 
evaluation of q iT (i* j s ). From the orthogonality of spher- 
ical harmonics we find that 



(0,0|17|3,r) = 



1 



47T 



(3.16) 



where r = 0, (to, c), (to, s), to = 1 — 3, and zero other- 
wise. The matrix quantities l|3.15|) were first introduced 
by Condon and Shortley for the description of atomic 
spectra^S but they are also at the center of the calcula- 
tion of the crystal electric field effects^ 



S,ja.4-f,ja 
I 

The general expression is 



In Eq. (|3~T4l> 



a,o— c.e 



stands for a radial average. 



dr r 2 J dr' r' K a (r)K b {r) 

K c (r')K e {r')vi{ry), (3.17) 
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where lZ a , T^-b, TZ-c and lZ e are radial components, and 
vi(r, r') — v\\(r, r') is given by Eq. I|3.1U[) . For the 7s5/ 3 
configuration the indices a, 6, c, e refer either to 7s or to 
5/ electron radial components. There are only two types 
of radial integrals, corresponding to nontrivial multipolar 
terms (I ^ 0) of 7s5/ 3 , which are 



v 



•v 



ff-ff 



drr 2 I dr' 'r'" 2 ' n 2 f (r)TZ 2 f (r') vi{r,r'),{^.18&) 



f»-sf _ 

V ~ 



(3.18b) 



drr 2 / dr'r' K s (r)K f (r) 



TZf(r')TZ s {r') vi(r, r'). 



Here I — 2,4,6 and V — 3, as follows from the selec- 
tion rules for c\(ff) and c\'(fs). The radial integrals 
v ff-ff anc j y/s-s/ are p r0 p 0r tional to the quantities F) 
and Gy introduced by Condon and Shortley in Ref. |2(j. 
It is important to notice that even the 7s electron with 
the trivial dependence of its angular part is strongly cou- 
pled to the three 5/ electrons through / — s and s — f 
transitions. 

The classification scheme for electronic transitions 
which we have introduced here is very useful for han- 
dling the single particle interactions as well. The main 
difference is that now the interaction occurs to a sin- 
gle electron while the rest of them produce Kronecker 
factors, Ref. This group of interactions includes the 
spin-orbit coupling H so , the crystal electric field Vcf and 
the mean field Vmf- The latter two interactions are dealt 
with in Sees. IV and V. The spin-orbit coupling is 



(3.19) 



where the sum runs over all electrons, V so being the 
corresponding one-electron spin-orbit operator. The 
s— electron does not experience the spin-orbit coupling 
and in the 7s5/ 3 case the summation includes only three 
5/ terms V so (i^), i* = 1 — 3, where 



dfl'(ft'), we obtain 

(I\n(I^V(R(n),R'(n'))\J') ri '\J)n = ^ 

4 4 

a(n)a'(n') a=l a'=l AA' 



13=1 



0'=1 



(3.21) 



Here the sum over a means the summation over all per- 
f f f 

mutations of indices Ji ,J2 ,,73 , J s at site n transform- 
ing them to indices j® (a = 1—4). Analogously, the 
sum over a' implies the summation over all permutations 

°f fx ji'Cji 3)i' S at site n' transforming them to j'* 
(a = 1—4). P(a) and P(a') stand for the parities of 
the permutations. Indices a and a' indicate which elec- 
tron at site ft interacts with which electron at site ft'. 
The other electrons labeled by (3 = 1 — 3 at site ft and 
by /3' = 1 — 3 at site ft' do not contribute to the in- 
teraction and produce the Kronecker delta symbols. The 
coefficients c\ are defined by Eq. H3.15fl . and the intersite 
(ft 7^ ft') interaction element \, a is given by 

vl a (ft - ft 1 ) = J dr r 2 J dr' r' 2 

xn 2 a (r)TZ 2 a ,(r')v AA ,(n,ftf; r,r'). (3.22) 

For the four electron space of 7s5/ 3 only even I and I' 
(in A, A') are retained in Eqs. (|3.21|) and Q3.22p . and 
= 0,2,4,6. Two very important examples of inter- 
site Coulomb interactions, namely, the crystal electric 
field and the mean field will be considered in sections IV 
and V. 



V so {if) = QL(it)S(it) 



(3.20) 



Here L(v ) and S(i^) are the one-electron operators of or- 
bital and spin momentum, Q is the constant of the spin- 
orbit coupling. The full intrasite Hamiltonian is given by 
Hintra = + H so . It describes the 7s5/ 3 configura- 
tion of a free neptunium ion. Since the present method is 
not based on perturbation theory it extends the classical 
calculations of Condon and Shortley 20 



E. Intersite matrix elements 

We start with expression (|3.8|l and write it in the space 
of many electron basis vectors |7), Eq. I|3.2II . Carry- 
ing out the angular integrations dfi(ri), dQ'(ft), dVt(ft') 1 



IV. CRYSTAL AND MEAN FIELD OF THE 5f 
NEPTUNIUM CONFIGURATION 



In this section we study the model where we assume 
that there are only three localized 5/ electrons at each 
neptunium site. Although we believe that the model is 
not adequate for a realistic description of Np02, espe- 
cially in the part concerning the loss of the magnetic 
moments in the ordered phase, it is nevertheless very in- 
structive to consider it in detail. The 5/ 3 configuration 
being relatively simple offers an opportunity to be stud- 
ied thoroughly and to understand the interplay between 
the disordered and ordered phases, or between the crys- 
tal and quadrupolar mean field. On the other hand, the 
configurations 7p5f 3 and 6c?5/ 3 involve too many basis 
states and consume to much time to be processed self- 
consistently for any temperature. 
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A. Free ion electron energy spectrum 



The basis states for the 5/ 3 configuration are given by 



\i 1 , i 2 , 



f\ 



(4.1) 



where as before, v = 1 — 14. The total number of ba- 
sis vectors is 14 x 13 x 12/3! = 364. We treat the 5/ 3 
configuration in the way which was specified in Sec. III. 

There are only / — / transitions described by four ra- 
dial integrals (|3.18a[) : v{ f ~ ff , I = 0,2,4,6. The others 
are zero due to the selection rules imposed by the co- 
efficients c A (ij), Eq. (|3.15l) . The radial integral Vq^~^ 
(Hubbard U) is not important here since it does not re- 
sult in term splittings. In the following we will use the 
condense notation F for // and thus = v F ~ F . 

These quantities are connected with the Slater (Condon- 
Shortley) parameters^ F'(5/, 5/) through the following 
relation: 



F-F 



47T 

21 + 1 



-F . 



(4.2) 



In particular, the Slater parameters F 2 , F 4 , and F 6 



of Amoretti et 



correspond to v F F = 14.007 eV, 



, F ~ F = 7.091 eV and v£~ F = 3.168 eV. (In order to 



obtain the exact term splitting quoted in Table IIVI of 
Ref. ^3 we nacl to scale their Slater parameters by a 
factor 0.9755.) Alternatively, the quantities v F ~ F can 
be calculated by using the radial dependence of the 5/ 
electrons IZf, Eq. l|3.18aj) . We have done such a calcu- 
lation and then corrected the parameters by comparing 
the splittings of the / 3 configuration with experimental 
data for Pr 3+ and Nd 4+ (details are given in Appendix 
B). We arrived at 



F-F 



18.164, 



F-F 



8.578, 



F-F 



Cf = 0.2547, in eV. 



3.362, 
(4.3) 



After calculating the matrix elements of the Coulomb 
repulsion and the spin-orbit coupling, we diagonalize the 
matrix 



Hi, 



y (3) + h,, 



(4.4) 



and obtain the electronic spectrum of 5/ 3 . The 9 low- 
est and two highest eigenvalues are shown in Table llVl 
where for comparison we also quote the spectrum of free 
Np ion of Ref. 43. Notice, that in comparison with the 
spectrum^ of Amoretti at al, 4 F 3 / 2 and 4 F 5 / 2 are higher 
than 4 /i3/2 and 2 Hg/ 2 , which is in better agreement with 
the sequence of terms for the 4/ 3 configurations of Pr 3+ 
and Nd 4+ known from experiment ^ 



B. CEF excitations in the disordered phase 

(T > 25 K) 

In the disordered phase (T > 25 K) the electron den- 
sity of the 5/ 3 configuration adopts the cubic (Oh) site 



TABLE IV: The 9 lowest and 2 highest eigenvalues of 5/ J , 
calculated with v F ~ F , Eq. 114. 311 . g stands for the Lande fac- 
tor. E A refers to the calculation of Amoretti et al, Ref. II2I 
Two highest values of E A marked by an asterisk were repro- 
duced by our calculation with the parameters of Ref. l42l 





term 


deg. 


9 (pb) 


E, meV 


E A , meV 


1 


4 r 


10 


0.7546 








2 


4 r 


12 


0.9704 


635.3 


657.7 


3 


4 r 

-'13/2 


14 


1.0993 


1204.2 


1256.2 


4 


F3/2 


4 


0.6027 


1244.2 


948.3 


5 


Hg/2 


10 


1.0154 


1617.7 


1438.3 


6 


F5/2 


6 


1.0067 


1702.9 


1399.3 


7 


4 T 

-'15/2 


16 


1.1797 


1715.3 


1762.3 


8 


£3/2 


4 


1.6546 


1861.6 


1614.9 


9 


-FV/2 


8 


1.1195 


1955.3 


1697.1 


40 


2 F 7/2 


8 


1.1317 


7796.6 


6541.1* 


41 


2 F 5/2 


6 


0.8589 


8008.2 


6631.6* 



symmetry. This density modulation is induced by the 
cubic crystal electric field experienced by three 5/ elec- 
trons. In terms of the multipole intersite expansion l|3.8|) 
it implies that for a given Np site n we treat 12 Np neigh- 
bors (ri[ = 1 — 12) and 8 oxygen neighbors (n' 2 = 1 — 8) in 
spherical approximation, i.e. I' = and Sa' (n') reduces 
to Yq = 1/VItt. The coefficients ca', Eq. Q3.15p . become 
simple, 



Co(ia<ja>) 



1 



S(ia',ja') 



(4.5) 



Here we write for A' = (I' = 0, Ai g ). At the central 
site n we expand the CEF in terms of SAF's S'a^H), 
Ai = (I, Ax g ), where Ai g stands for the unit representa- 
tion of the cubic site group Oh- The selection rules for 
the coefficients c\(i a ,j a ) of the / — / transitions imply 
that there remain only two nontrivial functions Sa ± with 
/ = 4 and / = 6, which correspond to the cubic harmonics 
K^(fl) and Kq(£1). The multipole two-center expansion 
(13.81) becomes 

V(R(n), R'(n')) = -±= V v Al0 (n, ft'; r, r 1 ) S Al (n), 

(4.6a) 



where 



v Al0 (n,n'; r,r') = 
1 



dn(n) dn'(n') 



S Al (n) 



\R(n)-R'(n')\ 



(4.6b) 



Here v Al o(n,n'; r,r') has a same value for all 12 Np 
neighbors (n^ = 1 — 12), and a same value for all 8 oxy- 
gen neighbors (n' 2 = 1 — 8). As follows from Eq. (|3.11() 
v Al o(n,ri' '; r, r') is independent of r' . Eq. I|3.22J) then can 
be written in the following form: 



Qfl a a 
V A! 



(n 



,,00 
v Al 



(4.7) 
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where 



and 



«a?o = / drr 4 Ki t (r)vAM*>**; f) ( 4 - 8a ) 



Qq' = J dr'r ,2 n 2 al (r'). 



(4.8b) 



Here the integrations are taken over < r' < Rmt, 
where Rmt is the radius of the muffin-tin sphere of nep- 
tunium or oxygen. (The influence of the interstitial re- 
gion will be discussed later.) Q a i refers to an electron 
at site ft' which interacts with one of the three 5/ elec- 
trons at n. We then can perform a summation over all 
electrons at fV and include also to this term the interac- 
tion with the nucleus. This results in replacing Q a i by 
zQmt in Eq. Ij4.7|l . Qmt and e being the total charge 
inside the MT sphere and the electron charge (e = —1), 
respectively. From Eq. i|3.11|l it follows that 



v a"o = «Ai a(n, ft'; Rmt, /') 



R 



MT 



where I in the index Ai is 4 or 6, and 



dr': 



• (1+2) 



(4.9a) 



(4.9b) 



Therefore, the CEF operator for any neptunium site 
(r < R M P T ) can be written explicitly as 



V CP (R(H))= Yl B i S (LA lg )(n) ( 

(=4,6 V 



R 



Np 
MT , 



where 



Bi = B 



Np 



B 



l > 



(4.10a) 



(4.10b) 



and 



B 



BP 



Np 



I 



i e f f evA 1 o{n, l n' l ;R M p T , /'), (4.10c) 



Q° f f ev Al0 (n, H' 2 ;R M P T , /'). (4.10d) 



We quote all relevant parameters of CEF in Table IV! As 
given by Eq. (|4.10a|l . the CEF operator Vcf is a one- 
electron quantityiMi 4 ^ CEF acts along with the Coulomb 
intrasite repulsion, Eq. I|4.4|l . Therefore, the total Hamil- 
tonian for the disordered phase becomes 



H dls (n) = H mtra + V CF (n). 



(4.11) 



Although we have considered CEF from first principles 
there is still an ambiguity related to the charge distribu- 
tion in the interstitial region. A more rigorous treatment 
of the problem is given in Refs. |4j47l A careful consid- 
eration of the problem based on the solution of a peri- 
odic Poisson's equation leads to a renormalization of the 



Np Np 
"l 







TABLE V: Calculated parameters of the CEF. 
fAi o(n, ft'i', -Rmt' ^ — is the central Np site, n[ is one 
of 12 Np nearest neighbors. v ( p q = Wa x o(n, n'2', R-mt' /'): 
n' 2 is one of 6 oxygen neighbors. 



v. 



Np o 
I 







Np Np 



B?/Q° f 

B Np ' 



5/p 



units 


/ = 4 


I = 6 


meV 


-816.7 


209.7 


meV 


-26.37 


-6.190 




0.1592 


0.0994 


K 


-3405.2 


546.0 


K 


-165.0 


-24.2 



charges inside the MT spheres, Ref. |47|. In other words, 
in Eqs. H4.1Uak O the effective charges for Np and O are 
given by 



ieff 



-AttR 



!MT — ■ 
MT 2_j 



0) 



3\{KR M t) 
K 



(4.12) 



where pj (K) is the Fourier series expansion of the elec- 
tron density in the interstitial region, is the spherical 
Bessel function. pi(K = 0) is the average density in the 
interstitial region, Q ou t/V out , where Q out and V out is the 
charge and volume of the interstitial region. 

The calculation of Q e f / according to Eq. 14.12|l is quite 
laborious since it requires the evaluation of the Fourier 
coefficients pj(K). Instead, below we consider two ap- 
proximation to (|4.12|l . In the first approximation we 
assume that Q^ f p f (I) = Q M P T and Q° ff (I) = Q MT , 

where Qmt anc ^ Qmt are ^ ne f°t a l charges inside the 
MT spheres of neptunium and oxygen. Here the electron 
charge in the interstices is completely ignored. In the 
second approximation we take 



eff(H) = Q 



MT 



0). (4.13) 



This expression corresponds to the homogeneous electron 
density distribution in the interstitial region. However, 

the modification of effective charges in this approxima- 

,n Pi 



and Q eff \ 



tion is too strong: Q e P j{I) = +4.108| 
+5.337|e| for neptunium, Q^^(I) 
Q?ff (II) = +0.945|e| for oxygen. As was discussed in 



^ Np f {II) = 
-0.284|e| and 

,(// 

the CEF splitting is overestimated in the sec- 
ond approximation. In reality the charge density in in- 
terstices is highly inhomogeneous, concentrated mainly 
in the proximity to oxygen and neptunium. This leads 
to pi(K ^ 0) ^ 0, and the last term in Eq. (|4.12|) acts 
in the opposite direction, decreasing Q e ff backward to 
Qmt values, which correspond to the first approxima- 
tion. 

The exact calculation of Q e f t according to Eq. I|4.12|l is 
beyond the scope of the present study. Instead, we have 
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FIG. 2: Splitting of the lowest 4 Iq/2 terms of the 5/ 3 con- 
figuration of Np in cubic crystal field, Eqs. l^.lOat d) as a 
function of the effective charges of neptunium and oxygen, 
Eq. J4.14al b1, Zero corresponds to the energy of the 4 Ig/2 
level of free Np ion. 

studied the crystal field effects as a function of Qeff by 
introducing 

Qeff( x *ff) = Qm P t + Xeff(Q? f P f (H) - Qmt)> (4- 14a) 
Q° ff (x eff ) = Qmt + x e ff(Q°ff(II) ~ QSt), (4.14b) 

where < x e ff < 1. Diagonalizing H dls , Eq. 14.11(1 . we 
have found that 41 terms of 5/ 3 are split into 120 distinct 
sublevels of Tq , I\ and r 8 symmetry of the cubic double 
group 0' h . In particular two lowest atomic- like levels are 
split according to the following scheme, 

% /2 ^r 8 + r 8 + r 6 , (4.i5a) 
4 / 11/2 ^r 8 ±r 7 ±r 6 ±r 8 . (4.i5b) 

The resulting splittings and the dependence of CEF on 
x e f f is shown in Fig. [21 The splittings of two lowest terms 
of the 5/ 3 configuration, Eq. l|4.15al bL is also given in 
Tables WT\ and IVTT1 for x e ff — and x e ff — 0.5, respec- 
tively. 

The most comprehensive study of the crystal field of 
the 5,f 3 configuration was performed by Amoretti et al, 
Ref . l42t Comparing their results with ours, we obtain the 
following relations connecting B4 and Bq with V4 and Ve 
used there, 

B 4 = (4.16a) 

B 6 = 16V8V 6 . (4.16b) 

We observe that for a realistic choice of Qeff, which cor- 
responds to x e ff ~ — 0.5, the CEF splitting is a few 
times smaller than the value 55 meV considered for CEF 
excitations in Ref. 0- Correspondingly, the calculated 



TABLE VI: CEF low energy spectrum and magnetic mo- 
ments of the 5/ 3 configuration of Np; x e jj = 0, Ae = 
7380.7 K. Calcula ted CE F parameters B 4 = -288.1 K, 
B 6 = 254.2 K, Eq. M.lOat . 





r 


deg. 


(U ~ ei) (K) 


Mz{[l B ) 




r 8 


4 





±(1.275, 1.429) 


4 r 
J9/2 


r 8 


4 


20.5 


±(0.517, 1.686) 




r 6 


2 


60.2 


±1.384 




r 6 


2 


Ae 


±1.778 


4 r 
J11/2 


r 8 


4 


Ae + 4.6 


±(1.241, 2.120) 




r 7 


2 


Ae + 6.3 


±1.775 




r 8 


4 


Ae + 26.1 


±(1.119, 3.835) 


TABLE VII: CEF low 


energy spectrum and magnetic mo- 


ments of the 5/ 


configuration of Np; x 


eff = 0.5, Ae = 


7396.3 K. 


Calculated CEF parameters B4 


= 1905.7 K, B 6 = 


-66.5 K. Ea. I|4.i0aj. 








r 


deg. 


(ft ~ ei) (K) 


Mz(hb) 




r 6 


2 





±1.383 


4 T 

J9/2 


r 8 


4 


56.8 


±(0.451, 2.077) 




r 8 


4 


164.9 


±(1.204, 2.333) 




r 8 


4 


Ae 


±(0.392, 3.862) 


4 r 

hi/2 


r 7 


2 


Ae + 28.6 


±1.777 


r 8 


4 


Ae + 110.8 


±(1.290, 2.845) 




r 6 


2 


Ae + 132.6 


±1.785 



parameters B 4 and Bq (Tables IvTI and |VTT|) are smaller. 
We discuss a possible explanation to this fact in the Con- 
clusions, Sec. VI. Notice that it is not possible to relate 
the feature at 55 meV with the 4 /n/2 splittings because 
it is situated at much higher energy ~ 650 meV, Table 
IIVI Most likely, the experimental excitations at 55 meV 
refer to the valence electrons delocalized on the Np-0 
bonds, while the lowest CEF excitations of Np lie at a 
smaller energy ~ 6 meV, Tables I VII IVIII 

Finally, we would like to notice that our CEF calcu- 
lation is based on first principles, and in that respect it 
differs from the othersf 4 ^ 4 ^^ which use fitting from ex- 
periment to extract the CEF parameters. There a re also 
other technical differences. In contrast to Refs. 1491501 
we do not assume that the full momentum J is a good 
quantum number, that allows for a mixing of components 
belonging to different J. In Ref. the basis was trun- 
cated to the first 11 low- lying levels. In our approach we 
do not have these limitations. 



C. Mean field and the structural phase transition 
at 25 K 

Now we consider the intersite quadrupole interactions 
VQQ(n,n') between a central Np site (sublattice {rti}) 
and its 12 nearest Np neighbors belonging to sublattices 
{n p r}, p' = 2, 3, 4. (Here we will not take into account the 
interactions involving higher spherical harmonics because 
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they are considerably smaller, Eq. i|3.11|l .') We will han- 
dle the Pn3m and Pa3 symmetries together since both 
of them result in the same local quadrupolar function 
(|2.1a|) for neptunium. 

The problem for the Pa3 structure has been considered 
m Refs.HHl where it is shown that this spatial order of 
quadrupoles gives an effective attraction between them. 
This direct quadrupole-quadrupole coupling can be cal- 
culated from first principles. However, we will see that its 
strength is not sufficient and we have to assume a sub- 
stantial reinforcement via oxygen mediated interaction. 
For the Pn3m structure the direct interaction is repul- 
sive and we have to resort to the indirect coupling from 
the very beginning. The important fact which we exploit 
in this section is that the rhombohedral (trigonal) mean 
field for both structures can be described by the same ex- 
pression (|4"2ljl below. 

We start by deriving an effective mean field for di- 
rect Coulomb coupling between quadrupoles in the Pa3 
structure. Following Ref.[lk, one obtains from Eq. (|3.21l) 
the following expression for the quadrupolar interaction 
operator between neptunium sites at ft\ and n p i 

^(n!,rv) = pfin,) pf(ft p ,). (4.17) 

Here fix € {«-i}, ft p > € {n p >}, and the quadrupolar den- 
sity operator at site n p is given by 

3 2 
I, J a a=l 13=1 

(4.18) 

f f f 

where again a is a permutation of j{ , j 2 1 3s transforming 
them to a new order given by j", n = 1 — 3. P{a) is 
the parity of the permutation, the index a stands for the 
interacting electron at the site. (The second interacting 
electron belongs to a neighboring Np site.) The other 
(noninteracting) electrons ((3 = 1,2) at n produce the 
product of the Kronecker delta symbols. The quadrupo- 
lar f — f coefficients arc defined as 

Cp(i a ja) = (ia\S p \j a ). (4-19) 

There are four types of such coefficients (i.e. p = 1 — 4) 
as follows from Eqs. (|2.1al -d'). Finally, the interaction 
constant 7^ ' in Eq. I|4.17[l is given by 

j ff = J dr r 2 J dr' r' 2 Tlj (r) TZ 2 f (r') v AA {ft, ft'; r,r'), 

(4.20) 

with v\\(jt,n'; r,r') defined by Eq. i|3.9f) . where ft — 
(0,0,0), ft' = (a/2)(0, 1,1) and A = (I = 2, T 2g ,k = 1). 
The corresponding SAF is Sa = Y 2 , Using the property 
(|3.11() we rewrite 7^ ' as 

// = _d v Q Q ( R Np i) _d_ (421) 



TABLE VIII: Calculated parameters of the mean field. 



v QQ {R N J T ) 


Q2/( r mt) 2 




X ff 


4568.2 K 


0.3095 


437.6 K 


1750.5 K 



where the short notation v ( ^ ( ^(R A ^ r ) stands for 
VAA(ft, ft 1 ', R-mti -^mt)i an d t ne "quadrupole charge" 
q{ =2 of 5/ electron is given by Eq. (|4.9bl) . 

In the mean-field approximation after summing over 
12 nearest neighbors belonging to 3 sublattices we arrive 
at the effective bilinear quadrupole-quadrupole operator 

U^(n p ) = -\ff (p])pj(n p ), (4.22) 

where \? f = 47^ > and the calculated parameters 
are quoted in Table IViTTl The same expression holds for 
PnSm, but in that case for the direct quadrupole inter- 
action < (repulsion), Appendix A. (p9) stands for 
an expectation value of the quadrupole operator Ij4.18|l . 
At zero temperature it is the quantum average 

(pf) = (I gs \pj(n P )\I gs ), (4.23) 

where \I gs ) refers to the ground state of the full mean- 
field Hamiltonian 

H MF (n) = U QQ (n) + V CF (n) + H mtra (n). (4.24) 

The intrasite part of the interactions Hi ntra is given 
by Eq. (|4.4() . For CEF we used the values x e ff = 

and Qeff — Qmt> Q?ff ~ Qmt- we discussed 

in Sec. IV. B this gives the most realistic estimate for 
CEF.) For both Pa3 and Pm3n structures the mean-field 
Hamiltonian has the Sq — C3 x i trigonal common site 
symmetry, and lifts the 4-fold degeneracy of the quartet 
states, 

T S ^E + E, (4.25) 

where E stands for the two-fold degenerate irreducible 
representation of C3. The states T 6 remain unsplit as a 
consequence of the Kramers theorem, r$ — > E. 

We now obtain a system of equations which can be 
solved self-consistently. At first we introduce an average 
(p9), which defines the interactions l|4.22|) . (|4.24() . In 
the space of the 5/ 3 configuration (364 state vectors) , we 
diagonalize the total Hamiltonian H MF , Eq. Il4.24[l . and 
obtain the eigenvectors \K), K = 1 — 364: 

H MF \K) =e K \K), (4.26) 

where the lowest value of ck corresponds to K = 1,2. 
This is the Kramers doublet of the ground state. We 
then calculate the quantities 

pf(K) = {K\pQ(ft p )\K}, (4.27) 

which evaluate the quadrupolar moments of the states 
K. Next, we find an improved value for (p?) which is 




1.0 1.1 



FIG. 3: A typical evolution of the order parameter amplitude 
(py } with temperature; po = {/9?}|t=o, Tq is the transition 
temperature. 
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FIG. 4: The transition temperature Tq as a function of the 
quadrupolar interaction constant A e / / . 



TABLE X: Mean-field (trigonal) splittings at T = with an 
enhanced quadrupolar interaction constant Ae/z/A^ = 10. 



TABLE 


IX: 


Mean-field (trigonal) splittings 


at T = 0, 




= 24.6 K, po = 


-0.0563. 








' = 1. 


T Q = 0.44 K, po = 


-0.0205. 














r 


deg. 


(ci-ex) (K) 


M z (pb) 


Q 
Pf 






(ci-d) (K) 


-M z (p B ) 




r 


deg. 


E 


2 





±1.9211 


-0.0563 




2 





±0.9464 


-0.0205 


E 


2 


54.0 


±1.1723 


-0.0014 




2 


1.1 


±0.7962 


0.0135 


E 


2 


72.8 


±1.3684 


-0.0004 




2 


20.7 


±0.6406 


-0.0083 


E 


2 


85.8 


±0.6211 


0.0249 




2 
2 


21.7 
60.8 


±0.8428 
±1.3840 


0.0140 
0.0011 


E 


2 


122.1 


±1.4893 


0.0256 













pj(K = 1) = pJ(K = 2), Eq. 11031) . The procedure 

continues until the input and output values for (pj ) con- 
verge. As a consequence of symmetry the expectation 
value (p?) is independent of the sublattice {n p } chosen 
for calculations, but the Hamiltonian and eigenvectors do 
depend on the choice. This is because the quadrupoles 
have different orientations for different sublattices. 

For nonzero temperature T the mean-field equation for 
{pj ) becomes 

K 

where 

Z = ^e- £ < (4.28b) 
K 

The results of the calculations are quoted in Tables Hxl 
H and in Figs. 3, 4, 5. With X f f =1750.5 K, Table lVIlTl 
which corresponds to the direct quadrupole-quadrupolc 
attraction of the Pai structure, we found that the transi- 
tion temperature Tq is only 0.44 K. A typical dependence 
of (pf) is shown in Fig. 3. The phase transition is of first 
order, with a discontinuity of the order parameter am- 
plitude (p?) = —0.0075 at Tq. Comparing the present 
calculation with that for cerium^ we observe that the 



low value of Tq is due to a small quadrupolar suscepti- 
bility of the ground state quartet Tg , since the difference 
in A is not that much (for cerium \ff = 2241 K.) 

The calculated transition temperature is very small in 
comparison with the experimental value of 25 K. We 
then conclude that in the framework of the 5/ 3 model 
the structural phase transition can not be explained by 
means of the direct bilinear quadrupole coupling. Most 
likely, there is indirect coupling between 5/ electron 
densities on neptunium via the oxygen atoms (superex- 
change interaction^). The microscopic consideration of 
the superexchange interaction is beyond the scope of the 
present work. Instead, we model it by increasing the 
value of \* f . In such a case A-^ becomes a phenomeno- 
logical parameter which we denote as A e //. By chang- 
ing X e f f we change the transition temperature as shown 
in Fig. 4. We have found that the experimental value 
of 25 K is achieved for \ e ffj\'> ~ 10, which indicates 
a substantial increase of the effective bilinear coupling, 
Eq. (|4.22|) . The relevant parameters of such strong mean- 
field are given in Table Ixl 

Finally, we would like to mention that the mean-field 
calculations have been done assuming that the CEF is 
weak, i.e. x e ff — 0. Increasing x e ff leads to an increase 
of the CEF splittings, which results in a strong suppres- 
sion of the transition temperature Tq, Fig.[S] Notice that 
at x e ff(r$ — > Tq) = 0.241 the ground state changes to 
the Tq doublet, Fig. [3 which is apparently unfavorable 
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-—•^ TABLE XI: The 5 lowest and the highest eigenvalue of 7s5/ 3 ; 

— - g is the Lande factor. 

* — : term deg. g (fiB) E 

=' : V 1 *h 9 06266 (T 



10 J \ 2 b h II 0.8983 196.5 



o- 

0.00 0.05 0.10 0.15 



\ : 3 5 h 13 1.0670 737.2 

X - 4 5 F 5 11 0.8700 814.4 

\ 5 5 I 7 15 1.1691 1245.1 

V 

125 30 82 1 F 7 1.0172 8229.5 



FIG. 5: The suppression of the transition temperature Tq 
with the increase of the CEF strength, x e ff, see also Fig. 2 
and Sec. IV.B; \ eff /\ ff = 10. 

for the quadrupolar order. With further increase of x e ff 
beyond the x e ff(Ts — > Tq) point the transition tempera- 
ture goes fast to zero and the ordered phase disappears. 



V. FOUR ELECTRON CONFIGURATIONS AT 
NEPTUNIUM SITE 

As follows from the previous section the model which 
takes into account only the 5/ 3 configuration at each nep- 
tunium site is not capable to explain the disappearance 
of the magnetic moments in the ordered phase. On the 
other hand, the charge distribution inside the MT sphere 
centered at the neptunium nucleus indicates that there 
is always approximately one valence electron instanta- 
neously present along with the three localized 5/ elec- 
trons, Table ITTT1 Even if the electron is in the s-orbital 
state, it experiences strong coupling with the 5/ electrons 
via intrasite s — f transitions. Therefore, the excitation 
spectrum at each neptunium site differs from that for 
5/ 3 considered in Sec. IV. In the following we model the 
couplings with the valence electron by considering 7s5/ 3 , 
7p5/ 3 , and 6d5/ 3 instantaneous configurations. Here we 
will not study the mean-field in such detail as for the 5/ 3 
configuration. Our main objective is to show that the 
ground state can be nonmagnetic and separated from the 
magnetic excitations by an energy gap larger than 25 K. 

A. 7s5/ 3 configuration 

The 7s5/ 3 configuration has been considered in detail 
in Sec. III. We are working here in the space of 728 state 
vectors \I) = |i{,z{,Zg; i s ). The 5/-states are coupled to 
two s— states through the f — s transitions accompanied 
by the multipole Coulomb interactions with I — 3. The 
strength of this interaction was estimated from the LDA 
calculation of a Np atom, Eq. (|3.18b|) . 

v( s ~ sf = 1.313 eV. (5.1) 



TABLE XII: CEF low energy spectrum and magnetic mo- 
ments of the 7s5/ 3 configuration of Np, Ae = 2293.2 K. CEF 



parameters B4 — 


-288.1 K, B a 


= 254.2 K; x eff 


= 0. 


r 


deg. (e 


i-ei) (K) 




E 


2 





0; 




3 


15.6 


±1.5671; 


Ti 


3 


25.1 


±0.2946; 


A 


1 


61.6 





Ti 


3 


Ae 


±2.2435; 


5 / 5 E 


2 


Ae+4.1 


0; 


T 2 


3 


Ae+6.5 


±2.2393; 


Ti 


3 


Ae+20.3 


±2.7113; 



The electron energy spectrum of the 7s5/ 3 configura- 
tion consists of 82 distinct levels, see Table IXll In the cu- 
bic crystal field two lowest levels and 5 is) are split as 
quoted in Table IxTTl Notice that the ground state is non- 
magnetic. However, at higher temperatures (T > 25 K) 
two low lying excitations of T2 and T± symmetry con- 
tribute to the Curie law for the magnetic susceptibility. 6 

In the ordered phase the local symmetry of the Np site 
changes to Sq (PaS) or D 3 d (Pn3m). As we discussed in 
Sec. IV. C the mean-field can be expanded in a multipole 
series and from symmetry it follows that the strongest 
interaction is of quadrupolar type, 

UQQ(n p ) = -Af p®(n p ), (5.2) 

where A^ = A // (p^) [compare with Eq. The 
quadrupolar operator p9(n p ) belonging to the sublattice 
{n p } is given again by Eq. (|4.18|l with the corresponding 
orientational function S p , Eqs H2.1al d). The only differ- 
ence is that the s— electron produces an additional Kro- 
necker symbol, i.e. in Eq. (|4.18() (3 = 1 — 3. Notice that 
both CEF, Eq. l)4.10al) . and the mean field, Eq. lf5T5)l. 
act only on the 5/ electrons. Taking A? — 5252 K we 
diagonalized the full Hamiltonian 

H MF (n) = U QQ (n) + V CF (n) + H mtra {n). (5.3) 

The resultant electronic spectrum is shown in Table IXTLTl 
We observe that the ground state is nonmagnetic, while 
the first magnetic excitation (E) does not contribute to 
the magnetic susceptibility if T < 25 K. 
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TABLE XIII: Mean-field (trigonal) splittings of 7s5/ 3 at T = 
0, A f = 5252 K. 



r 


deg. 




-a) (K) 


Mz(hb) 


A 


1 










E 


2 




26.0 


±0.3467 


E 


2 




124.1 


±0.6786 


A 


1 




301.8 





A 


1 




330.1 





E 


2 




516.5 


±1.4769 


TABLE 


XIV: The 5 lowest and the highest eigenvalue of 


7p5/ 3 ; g 


is the Lande factor. 








term 


deg. 


9 Q*b) 


E 


1 




11 


0.7430 





2 


5 h 


9 


0.7459 


154.5 


3 


5 K 6 


13 


0.9450 


689.7 


4 


5 h 


11 


0.9875 


756.2 


5 


5 G 2 


5 


0.5948 


1248.1 


242 


3 Pi 


3 


0.5031 


10515.7 




B. 


7p5/ 3 configuration 





TABLE XV: CEF low energy spectrum and magnetic mo- 
ments of the 7p5/ 3 configuration of Np, Ae = 1789.1 K. CEF 



parameters B4, — — 


288.1 K, B 6 


= 254.2 K; x eff 


= 0. 


r 


deg. ( 


u - ei) (K) 


Mz{^b) 


T 2 


3 





±1.8525; 


5 K 5 Ti 


3 


3.2 


±1.7788; 


E 


2 


16.7 


0; 


Ti 


3 


45.0 


±2.1534; 


E 


2 


Ae 


0; 


5 h T 2 


3 


Ae±14.8 


±1.8645; 


Ti 


3 


Ae±24.8 


±0.3655; 


A 


1 


Ae±59.1 





TABLE XVI: Mean-field (trigonal) splittings of 7p5/ 3 at T = 


0; A f = 6612 K, A p 


= 3426 K. 






r deg. 




- ft) (K) 


M z {^b) 


A 1 










E 2 




24.6 


±0.4171 


E 2 




105.8 


±0.8134 


A 1 




246.4 





A 1 




272.6 





E 2 




480.7 


±1.7473 



In case of the 7p5/ 3 configuration, we construct (14 x 
13 x 12/3!) x 6=2184 basis vectors 

\I) = \i{,44;i"). (5.4) 

As before, v refers to a 5/ electron, i p = (k, s z ) to the 
7p electron (fc is its orbital part, k = 1 — 3, and s z is the 
spin part). 

The interaction between 5/ electrons was described in 
detail earlier. In addition, p — p and / — / transitions 
lead to intrasite multipole interactions with the I = 
and 1 = 2 components, f—p transitions give the Coulomb 
multipole couplings with I = 2 and I = 4. The relevant 
Slater parameters were extracted from the radial depen- 
dences of a Np ion, 

v[ v = 1AA2, v{ p = 0.739, 

v^-pp = 5.237 (in eV). (5.5) 

Also, the p-electron experiences the spin-orbit interaction 
with C P = 1.2795 eV. 

First, we calculated the electron spectrum of the free 
ion, Hi n tra and obtained 242 distinct levels, Table IxTVI 
In the cubic crystal the 5 K4 and 5 l4 levels are split as 
quoted in Table IXVI (The procedure of treating CEF 
effects is outlined in Sec. IV. B, CEF does not act on 
the p— electron.) Notice, that the ground state now is a 
magnetic triplet of the Ti symmetry, which together with 
two other magnetic T\ excitations, Table fXVl gives the 
Curie law for the magnetic susceptibility. 

At T c = 25 K the structural phase transition occurs 
and the symmetry of the Np sites is reduced. This sym- 
metry change is accompanied by lifting degeneracies of 



some cubic levels. In particular, the ground state Ti 
triplet is split in a doublet and a single level as demon- 
strated in Table IXVll The mean-field is approximated 
by its quadrupolar electric part, 

U^{fU) = -A^(^)-A*>pQ(^), (5.6) 

where 

A? = \ff(pf) + \fv(pQ), (5.7a) 

AP = \Pf(pf)+\PP(pQ). (5.7b) 

The following parameters of the interaction were as- 
sumed, 

p 

— P — = 0.1604, A f = 6612 K, A p = 3426 K. 

(5.8) 

Since the charge density expansion of the p-electron has 
quadrupolar components, Cf.(i p j p ) = (i p \Sk\j p ) ^ (Sk 
are given by Eqs l|2.1al -d)). the mean- field expansion 15. 6|) 
includes the quadrupolar projection on p-states, 

3 

I, J a /3=1 

(5.9) 

f f f f 
Here a is a permutation of j{ , J3 to j c 0, p = 1 — 3; 

P{a) is the parity of the permutation. The quadrupolar 

operator for the 5/ electrons is given again by Eq. (|4.18|) 
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TABLE XVII: The 5 lowest and the highest eigenvalue of 
6d5/ , g is the Lande factor. 





term 


deg. 


9 (Ms) 


E 


1 


b L 6 


13 


0.7530 





2 


5 K 5 


11 


0.7201 


167.2 


3 


B L 7 


15 


0.9232 


667.9 


4 


5 K 6 


13 


0.9269 


767.2 


5 


3 D 3 


7 


0.6800 


886.9 


383 


3 Pi 


3 


1.0004 


13116.4 



where the index (3 comprises the additional p— electron, 
i.e. 0= 1-3. 

Notice that in the ordered phase the ground state level 
is single and nonmagnetic, Table IXVII This mechanism 
can explain the loss of magnetic moments because the 
magnetic excitations of the 7p5f 3 configuration lie too 
high in energy. 



C. 6d5/ 3 configuration 

The basis vectors here are 

\I) = \i{A,i{;i d ), (5-10) 

where index stands for 5/ states {¥ = 1 — 14), while 
the index i d = {k, s z ) refers to five d-orbitals and the spin 
projection s z . Thus, i d = 1 — 10, and in total there are 
3640 nonequivalent basis vectors \I). 

We start by considering the intrasite interactions 
Hi n t ra . Here, in addition to / — / interactions we distin- 
guish two groups. The first group arises between d — d 
and / — / transitions. It is described by the multipole 
Coulomb repulsion with the I = 0, 2 and 4 angular com- 
ponents (SAF's). The second group is due to the / — d 
and d— f transitions. The corresponding multipole inter- 
actions are with I — 1, 3 and 5. The relevant parameters 
were extracted from the LDA calculation of the Np ion 
in the 6d5f 3 configuration, 

v fd-df _ n 322i v f d - d f = 3J01, v( d ~ df = 1.794, 

v ff-dd = n 2g9) v ff-dd = 3 4g2 ( m e y) ; 

Cd = 0.3497 eV. (5.11) 

The parameters for the f — f interactions were kept un- 
changed. We then diagonalized the 3640 x 3640 matrix 
of (I \H. i n tra\ J) and obtained 383 distinct levels. The 5 
lowest and the highest levels are quoted in Table IXVIII 
The CEF splittings of the lowest 5 Lq and 5 K$ levels are 
given m Table IXVIIll It should be noted that unlike be- 
fore, the CEF operator acts not only on the 5/ electrons 
but also on the Qd one^ 

U CEF (n) = B{ pj(rl) + B{ p 6 f (n) + B d p d (n). (5.12) 



TABLE XVIII: CEF low energy spectrum and magnetic mo- 
ments of the 6d5/ 3 configuration of Np, Ae = 1949.5 K. 



CEF 


parameters 


Bi = 


-288.1 K, B 6 = 


254.2 K, and 


Bf = 


-232.7 K; x 


eff = 0. 








r 


deg. 


(d - ei) (K) 


M z (pb) 




A 


1 










T 2 


3 


9.7 


±0.1809; 


5 r 


A 


1 


12.6 







T-l 


3 


15.6 


±0.3772; 




E 


2 


36.4 


0; 




T 2 


3 


43.6 


±2.0637; 




Ti 


3 


Ae 


±1.7253; 




E 


2 


Ae+3.1 


0; 




T> 


3 


Ae+21.1 


±1.8006; 




Ti 


3 


Ae+25.3 


±2.0883; 



Here p l f(n), 1 — 4,6, and p\{n) are cubic projectors on 
the / and d states, respectively, given by 

3 3 
/,./ a a=l /3=1 

(5.13a) 

p 4 d (n) = e \i) j2 p (^ d \ K ^ ad ) f[ S i^( J \' 

I, J a 0=1 

(5.13b) 

where Ki (O) refers to the cubic harmonics with I = 4 and 
6. Here we keep the same notations as before [Eqs (|4.18|l 
and l|5.9l) ]. i.e. the permutation a transforms the indices 

jiiji^ji t° j a t, K = 1 — 3. The permutations which 
interchange the d and / indices are excluded because they 
give zero contribution to l|5.13al bL The parameter Bf 
was calculated by the method described in Sec. IV. B. 
For x e f f = we have found that 

— ^ = 0.1287, Bf = —232.7 K. (5.14) 

\n MT ) 

Notice that for the ground state the CEF gives a non- 
magnetic single level, Table IXVlIII but there are two low 
lying magnetic levels (T2 and Ti) at 9.7 and 15.6 K, which 
contribute to the Curie law of the magnetic susceptibility 
at T > 25 K. 

Below 25 K the local symmetry of Np is lowered. The 
mean field is given by 

UQQ(n p ) = -Afpf(fi p )-A d pQ(n p ), (5.15) 

where 

tJ =\"{p ( })+\t*(p% (5.16a) 

K d = \ d f{p ( j)+\ dd {p ( j). (5.16b) 

Here again p?(n) and p^(n) are quadrupolar projection 
on the / and d states, respectively. They are given by 
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TABLE XIX: Mean-field (trigonal) splittings of 6d5/ 3 at T = 
0; A / = 6220 K, A d = 3442 K. 



r 


deg. 


(K) 


Mz{hb) 


A 


1 








E 


2 


32.8 


±0.4362 


E 


2 


136.6 


±0.8797 


A 


1 


311.2 





A 


1 


337.0 





E 


2 


601.2 


±1.7673 



expressions similar to <|5.13al -b). where we replace Ki(n p ) 
by S p , Eqs l|2.1al d). for 4 sublattices {n p } = 1 — 4 of Pa3 
or Pn3m, Sec. II. [Compare also with Eqs l|4.18|l and 
(15.91) .] Below we approximated the parameters of this 
interaction by 

— ^ — =0.1713, A / = 6220 K, A rf = 3442 K. 
« P t) 2 

(5.17) 

We then diagonalized the whole Hamiltonian H MF 
(H MF = U QQ (ft) + V CF + H intra ) and obtained the low- 
est energy levels quoted in Table IXIXl Notice that now 
the first magnetic excitation of E symmetry is separated 
from the nonmagnetic ground state by an energy gap of 
~ 33 K, which implies again disappearance of the Curie 
law for the magnetic susceptibility of the ordered phase 
(T < 25 K). 



VI. DISCUSSION AND CONCLUSIONS 

In resonant X-ray scattering (RXS) experiments the 
transition to the ordered phase at T c = 25.5 K in 
Np02 manifests itself by the appearance of superstruc- 
ture Bragg reflections [like (003) and others], which are 
not compatible with the FmZm structure of the disor- 
dered phase>ii*i£ Paixao et at, Ref- ITEl ascribed the sym- 
metry of the ordered phase to the Pn3m space group 
(No. 224, Ref. l22|L but did not elaborate on the theory. 
In this work (Sec. II), besides Pn3m we discuss also the 
space group Pa3 (No. 205), which is the other possible 
candidate for the triple-ff" quadrupole ordering. Pn3m 
and Pai are close symmetries. They imply condensations 
of different modes (X± and X%, Ref. l23f) at the same 
X— point of the Brillouin zone. We have considered the 
dependence of the scattering amplitude for different po- 
larizations on the azimuthal angle ip and the Bragg angle 
taking into account the domain pattern of both sym- 
metries. Pn3m and PaS produce superstructure Bragg 
reflections at the same sites (h,k,l) of the reciprocal lat- 
tice (Tables H] and |HJ . The relation between the space 
groups is such that in RXS experiments the (00£) reflec- 
tion of Pn3m behaves like the (£00) and ((W0) reflections 
of Pa3 (I = In + 1). Therefore, we believe that a spe- 
cial care should be taken to distinguish between these 



two symmetries. We have shown that the two structures 
are completely different with respect to the direct bi- 
linear quadrupole-quadrupole interactions. The Pn3m 
symmetry leads to a repulsion between the quadrupoles 
while the PaZ structure implies an attraction. (This is 
the main reason why the Pa3 symmetry occurs in many 
molecular solids: NaO^, Naffi, C 60l 26 . 28 i 29 . 30 i 31 ) 

We have presented an ah initio approach to crystal- 
and mean field for the structural phase transition at 25 K. 
The method is based on the technique of expanding the 
Coulomb repulsion between electrons in a multipolar se- 
ries, Sec. III. 

In the disordered phase we considered the Hamiltonian 
which includes the crystal electric field (CEF) effects and 
the intrasite Coulomb repulsion responsible for Hund's 
rules on equal footing, Sec. IV. The crystal electric field 
levels above the ground state are in fact the lowest ex- 
citations of the neptunium electron complex. The typ- 
ical splittings there are of the order of ten kelvins and 
thus the electron spectrum is very sensitive to the crys- 
tal symmetry. In this paper we have pursued the simple 
idea that the symmetry lowering produces a splitting of 
the many electron spectrum, which can explain the dif- 
ference in the behavior of the magnetic susceptibility and 
the loss of the magnetic moments. Our band structure 
calculation (Sec. III. A) indicates that besides the three 
localized 5/ electrons there is always approximately one 
conduction electron at each Np site. This changes the 
effective instantaneous configuration from the three elec- 
tron 5/ 3 to a four electron one (7s5/ 3 , 7p5/ 3 or 6c?5/ 3 ) 
and opens a possibility to obtain a nonmagnetic ground 
state without invoking the concept of the octupole or- 
der parameter £ From this point of view it represents an 
alternative to the latter, and we believe that both ap- 
proaches deserves a thorough experimental consideration 
and verification. 

Our finding that the CEF splitting in Np02 is rather 
small (~ 50 — 150 K) contrasts the commonly accepted 
value of 55 meVi 4 ^ We think that the excitation at 
55 meV observed experimentally by neutron scattering 42 
is a higher excitation of the neptunium electronic com- 
plex possibly involving the valence electrons on the Np-0 
bonds. The scale of CEF excitations calculated in the 
present work is a better match for the width of peaks in 
magnetic neutron scattering cross sections^ We believe 
that the crystal and mean field should be of the same or- 
der of magnitude and the transition temperature of 25 K 
gives a natural estimation for it. An implicit support of 
our viewpoint is the fact that in the ordered phase there 
appears an inelastic peak centered at about 6.4 meVi^ 

The phase transition to the ordered phase sets in at 25 
K and reduces the basic local symmetry of the neptunium 
sites to Sq = C3 x i (or 3). This symmetry holds for both 
PnZm and Pa3 space groups. [The addition of three two- 
fold axes or three mirror planes for Pn3m changes the 
point group to D 3 d (or 32) but this is not of principal im- 
portance.] We have found that the structural change can 
not be ascribed to the bilinear quadrupole-quadrupole 
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Coulomb repulsion, which is too weak to drive the tran- 
sition. Therefore, one has to resort to the superexchangc 
interactions^ via oxygen as a driving force for the trans- 
formation. In this work we have introduced an effective 
(enhanced) quadrupole interaction and studied the inter- 
play between it and the crystal electric field. The influ- 
ence of both interactions on the transition temperature 
has been investigated in detail for the 5/ 3 configuration, 
Sec. IV and figures 3, 4, and 5. 

The most intriguing question is the disappearance of 
the neptunium magnetic moments below 25 K. This how- 
ever can be explained if we consider a four electron com- 
plex at the neptunium site. We have demonstrated that 
all relevant configurations (7s5/ 3 , 7p5/ 3 , 6d5/ 3 ) can lead 
to a nonmagnetic ground state separated from the mag- 
netic excitations by an energy gap larger than ~ 25 K, 
Sec. V. Perhaps, the most clear example is the 7p5/ 3 con- 
figuration fTables lXVI IXVl(l . In the disordered phase the 
ground state is a triplet ( Tabic IXVl) . while in the ordered 
phase it becomes a singlet fTablc lXVIf) . The general idea 
for the loss of magnetic moments is similar to the one 
used by Kondo and Anderson and often referred to as the 
Kondo effectiSi Notice however, that here we are dealing 
with the intrasite interactions treated on ab initio level. 
In particular, we replace the Anderson hybridization 52 
which is linear in terms of creation/annihilation opera- 
tors for valence and localized electrons, by the Coulomb 
intrasite repulsion, which being a density-density cou- 
pling is bilinear in terms of these operators. Another im- 
portant theoretical ingredient of our model is the sym- 
metry lowering which modifies the excitation spectrum 
of the electron system at low temperatures. This part is 
absent in the Kondo mechanism. 

CEF and mean field have been objects of theoreti- 
cal interest for many yeara^^^SiS and we would like 
to mention here some important relations between our 
model and other approaches. We have shown that CEF 
effects can be perceived as a first meaningful term of 
the intersite multipole expansion, when all neighbors of 
a neptunium site are considered in the spherical approx- 
imation (I' = 0). It is then reduced to a single parti- 
cle potential— ii 5 . The intersite nonspherical terms are 
also included in the full potential (FP) electron band 
structure calculations like FP-LMTO (linear muffin-tin 
orbital method) and FP-LAPW (linear augmented plane 
wave method). 3 " Therefore, in principle one can say 
that the CEF effects are equivalent to the full potential 
treatment— iii However, there are two very important 
caveats here. First, in the band structure calculations 
the nonspherical terms of the potential apply to itinerant 
electrons in the ground state, while CEF effects are con- 
sidered usually for localized electrons in the ground and 
excited states. The second more important objection is 
that practically all band structure calculations are based 
on the single determinant approximation. This intrinsic 
feature does not allow to describe the intrasite interac- 
tions fully. In particular, the atomic term structure and 
Hund's rules are excluded from the consideration. This 



shortcoming does not apply to our treatment which is 
based on a many determinant approach. For the intra- 
site part of interactions our model is very close to the 
scheme described by Condon and Shortley for the elec- 
tron spectra of atoms and ions although there are some 
unimportant differences. The main approximation and 
limitation of our approach is the choice of basis. Once the 
many electron basis states are fixed, the work of calculat- 
ing matrix elements is reduced mainly to classifications 
of all possible permutations, which is done numerically. 
Thus, the fundamental group of electron permutations 54 
is explicitly taken into account. No additional approxi- 
mations used sometimes for crystal field calculations are 
made. This distinguishes our approach from Stevens— 
where the CEF is expressed in terms of equivalent oper- 
ators J x , J y , and J z . The latter approach as well as the 
work of Lea, Leask and Wolf for cubic CEF— based on it 
are justified only if J is a good quantum number. Notice 
also that the approach of Stephens starts with the sym- 
metry arguments while the interactions are introduced 
later in a phenomenological manner. 

However, the present calculation scheme does not take 
into account chemical bonding in an intrinsic way. There- 
fore, further development of the method should be fo- 
cused on this problem. 
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APPENDIX A: BILINEAR QUADRUPOLE 
INTERACTIONS FOR Pn3m AND PaS 

Here we calculate the direct bilinear electronic 
quadrupole-quadrupole interactions for Pn2>m and Pa3 
structures, Fig. We consider the quadrupolar compo- 
nents Si of T2 g symmetry (i = 1 — 3) at site ft = = 
(0,0,0) (ft £ {"-i})- There are 12 nearest neighbors of 
ft belonging to three different sublattices {n p }, = 2 — 4. 
The interactions between three components Si centered 
at ft = and those (Sj) located at four nearest neighbors 
(ft 1 = 1 — 4) of the second sublattice (ft 1 £ {1-2}) are given 
in Table IXXl 

The matrix S 2 '(n) — Sj(ft') for the fourth sublattice 
{ft.4} is given by the same Table provided that the cyclic 
permutation S\ — > S' 3 , S2 — » S[, S3 — > S' 2 is performed. 
Here we label four nearest neighbors ft' = 5—8 of ft as 5 = 
o(l/2, 1/2, 0), 6 = o(-l/2, 1/2, 0), 7 = a(-l/2, -1/2, 0), 
and 8 = a(l/2, -1/2,0). The matrix S'?(n) - S'j(ft') for 
the third sublattice {n.3} is obtained from Table IXXl bv 
replacing Si -> S%, S 2 -> S3, S 3 -> S". Here ft' runs 
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TABLE XX: The matrix of interaction Si(n) — Sj(h") between 
three quadrupolar components of T2 9 symmetry, 5*i = Y2 S , 



S 2 



S3 



Si are centered at n = (0,0,0), Sj at 



four neighbors n of the second sublattice {712}- 



n Coord. 



1 o(0, 

2 a(0, 

3 a(0, 

4 o(0. 



1 

■ 2 ■ 



1 i) 

? -11 
2' 2^ 

-I) 



(y) = (1,1) (2,2) (3,3) (1,2) (1,3) (2,3) 

7 a a /3 

7 a a — P 

7 a a /3 

7 q a — P 



For the Pa3 structure we have the function Sb, 
Eq. pTb)l . at n' = 1 - 4, and £^ Q ({n 2 }) = -4 7 /3. The 
same value is obtained for E^ Q ({n 3 }) and E^ Q ({ni}). 
As a result we arrive at 



E QQ (Pa3) = -4 7 < 0. 



(A3) 



We conclude that the Pn3m structure always leads to 
a repulsion while the Pa3 to an attraction. 



over the sites 9 = a(l/2, 0, 1/2), 10 = a(-l/2, 0, 1/2), 
11 = o(-l/2,0 L -l/2), and 12 = a(l/2, 0, -1/2). 

For the Pn3m order of quadrupoles we have iS a , 
Eq. fl2~Ta|) . at (0, 0, 0) and S d , Eq. fTTd|) . for n' = 1 - 4. 
By means of Table IXXI we find that this interaction is 
given by 

^o W (W) = ^(7-2a). (Al) 

The same expression is obtained for the interactions with 
four neighbors of third and forth sublattices, ({n 3 }) 
and £?q ({rc.4}). The total interaction then is 

EQQ(Pn3m) = E$ Q ({n 2 }) + E® Q \{n 3 }) + E$ Q ({n 4 }) 
= 4(7- 2a) > 0, (A2) 

because 7 > and a < independently of the lattice 
constant a. 



APPENDIX B: CORRECTION OF SLATER 
INTEGRALS FOR NP 

Since the experimental data on the energy splittings 
of the 4/ 3 configuration of Pr 3+ and Nd 4+ are available 
from the atomic database of NIST, Ref. l43l while there 
is no such information for 5/ 3 , we have performed cal- 
culations of v F ~ F (I = 2,4,6), Eq. (|3.18a|l . by using the 
radial dependence of IZf obtained from LDA calculations 
of atoms. After this we diagonalize the Hamiltonian of 
the free ion (V^ + H so ) and compared our calculated 
spectra with the experimental ones. We have noticed 
that the comparison is improved (the sequence of terms 
corresponds to the experimental one) if we reduce v F ~ 
and Vq~ f by a factor of 0.75 while keep v F ~ F almost the 
same (factor 0.975). Therefore, we have used the same 
scale factors for Np in Np0 2 and obtained parameters 
given by Eq. I|4.3|l . 
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